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Abstract 

We carry out a systematic investigation of all the 2-loop integrals occurring 
in the electron vertex in QED in the continuous D-dimensional regularization 
scheme, for on-shell electrons, momentum transfer t = —Q 2 and finite squared 
electron mass m 2 = a. We identify all the Master Integrals (Mi's) of the 
problem and write the differential equations in Q 2 which they satisfy. The 
equations are expanded in powers of e = (4 — D)/2 and solved by the Euler's 
method of the variation of the constants. As a result, we obtain the coefficients 
of the Laurent expansion in e of the Mi's up to zeroth order expressed in close 
analytic form in terms of Harmonic Poly logarithms. 
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1 Introduction 



The QED electron form factors at two loops were considered in pQ, for massive on- 
shell electrons and arbitrary momentum transfer within the Pauli-Villars regulariza- 
tion scheme and giving a fictitious small mass A to the photon for the parametrization 
of infrared divergences. The main results of [I] are the analytic calculation of the 
imaginary parts of the form factors for arbitrary momentum transfer in terms of 
Nielsen's polylogarithms 0111], and of the charge slope of the electron at two loops 
(besides the check of the magnetic anomaly). The analytic evaluation of the real 
parts, expected to involve a class of functions wider than Nielsen's polylogarithms, 
was not attempted in [T]. 

To our knowledge, the full analytic calculation of the real parts of the 2-loop 
QED form factors, for arbitrary momentum transfer and finite electron mass, has 
not yet been carried out, despite a great number of papers dealing with a variety of 
kinematical configurations (neglecting typically the electron mass at large momen- 
tum transfer). 

In this paper we work out a systematic investigation of all the 2-loop integrals 
occurring in the electron vertex in QED in the continuous Z)-dimensional regular- 
ization scheme j3] (using the same D for ultraviolet and infrared regularization) for 
on shell electrons of finite squared mass m 2 = a and arbitrary momentum transfer 
t = —Q 2 . We identify all the Master Integrals (Mi's) occurring in all the graphs and 
evaluate them analytically, in terms of Harmonic Polylogarithms [HUH], for arbitrary 
value of Q 2 . We present the results for spacelike momentum transfer, i.e. t < or 
Q 2 > 0; the case of timelike t can be obtained by standard analytic continuation. 

The extraction of the form factors and their expression in terms of the Mi's will 
be carried out in a subsequent paper. 

The diagrams involved are those shown in Fig. (jl]). 

Following a by now standard approach, we first express all the scalar integrals 
associated to each graph in terms of the Master Integrals (Mi's) by using the inte- 
gration by parts [7j and Lorentz invariance jH] identities, then write the differential 
equations on the momentum transfer which are satisfied by the Mi's, and expand 
the equations in powers of e = (4 — D)/2 around e = (D = 4) up to the required 
order. We obtain in that way a system of chained differential equations for the coef- 
ficients of the e-expansion of the Mi's and finally solve the system for the coefficients 
by Euler's variation of constants method. As a result, we express the coefficients in 
close analytic form in terms of Harmonic Polylogarithms j3J |U] . 

Let us recall that the Euler's method requires the solution of the associated 
homogeneous equation. Even if general algorithms for the solution of differential 
equations are not available, it is to be stressed here that all the homogeneous equa- 
tions which we had to solve came out to be essentially trivial (typically, first order 
homogeneous equations with rational coefficients). 

The present paper is structured as follows. In section El we review briefly the 
techniques for reducing the calculation of generic multi-loop Feynman graph inte- 
grals to the calculation of the Mi's. Integration by Parts, Lorentz Invariance and 
general symmetry relations are recalled and the application of this approach to our 
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Figure 1: 2- loop vertex diagrams for the QED form factor. The fermionic external 
lines are on the mass-shell p\ = p\ = —a, while the wavy line on the r.h.s. has 
momentum Q = pi + P2, with Q 2 = — s. 
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case is discussed. In section |3] we review the method of differential equations for 
the calculation of the Mi's. In section 0] we describe exhaustively the case of three 
typical integrals, giving in some details the system of differential equations and the 
steps for obtaining the solution. In section El we present the results for all the Mi's 
encountered in the calculation of the 2-loop vertex diagrams and in section |U] we 
give the results for the scalar 6-denominator vertex diagrams which are not Mi's. 
Sections [7| and |H1 contain respectively the expansions of the vertex 6-denominator 
diagrams in the region of great and low momentum transfer. Finally, after the sum- 
mary, section El and appendix El where we give the routing used for the explicit 
calculations, in appendix El we give the results of the 1-loop diagrams involved in 
our calculations and in appendix Owe list the results for all the reducible diagrams 
appearing in the calculation. 



2 The reduction to master integrals 

The aim of this paper is the evaluation of all the possible scalar integrals which can 
occurr in the calculation of the Feynman diagrams of Fig. ((TJ. They imply two loop 
momenta, k\ and k 2 , and three external momenta, Pi,p%,Q', among them only two 
are independent, because of the momentum conservation law: p\ + p 2 = Q- With 
two external momenta and two loop momenta, we can construct three Mandelstam 
invariant variables, p\-,p\ and Q 2 = (pi +P2) 2 (we use the Euclidean metric, so that 
the mass shell conditions are p\ = p\ = —a, where a is the squared electron mass) 
and seven different scalar products involving the loop momenta, namely (pi ■ k\), 
(Pi ■ h), (P2 ■ h), (P2 ■ h), (h ■ k 2 ), (kj) and (Af). 

The graphs of Fig. (£[]) involve up to 6 different propagators: more exactly, the 
graphs (a,b,c) have 6 different propagators, while graph (d) contains 2 equal electron 
propagators, graph (e) 2 equal photon propagators, so that two graphs, (d,e), involve 
only 5 different propagators. In the following we will not consider anymore graphs 
but topologies: topologies will be drawn exactly as Feynman graphs, except that 
all propagators are different. To make an example, when a graph contains twice 
some propagator, as the two equal photon propagators of graph (e) above, the 
corresponding topology contains that propagator only once; indeed, the topology 
of graph (e) of Fig. (JTJ) is given by the topology (e) of Fig. (jU). Besides those 
topologies, we will also encounter all the subtopologies obtained by removing from 
the graphs one or more propagators in all possible ways. 

Let t be the number of the propagators in any of the topologies or subtopologies; 
we can express t of the 7 scalar products containing the loop momenta in terms of 
the propagators, (the remaining (7 — t) scalar products will be called irreducible) 
and correspondingly the most general scalar integral associated to that topology or 
subtopology has the form 

{d D h}{d D k 2 }^-^, (1) 

where {d D k} is the loop integration measure (its explicit expression, irrelevant here, 
will be given in sectional), the integer m ; , i = 1, t are the powers of the t propagators, 
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with mi > 1, and the integer rij,j — 1, q, q = (7 — t), with rij > 0, are the powers of 
the irreducible scalar products. Let us further recall that the continuous dimensional 
regularization makes the definition meaningful for any values of the 7 integer rii, rrij. 

We will denote with J t rjS the family of the integrals with a same set of t propa- 
gators, a total of r = X]i( m « — 1) powers of the t propagators and s = ^ ■ rij powers 
of the (7 — t) irreducible scalar products. The number of the integrals contained in 
the family is 

"^(tt 1 )^: 6 )- (» 

As we will see more in detail in a moment, one can establish several identities 
involving integrals of the type of Eq. (JTJ) with different sets of the 7 indices mi,rij. 
The identities can be written in the form of a sum of a finite number of terms set 
equal to zero, where each term is a polynomial (of finite order and with integer 
coefficients in the variable D, a and the Mandelstam invariants) times an integral of 
the family, as will be seen explicitely in the example of next section. 

The identities can be used to express as many as possible integrals of a given 
family in terms of as few as possible suitably choosen integrals of that family - called 
the Master Integrals of that family. 

The identities will be generated by using Integration by Parts, Lorentz Invariance 
(or rotational invariance in D dimensions) and symmetry considerations. 

2.0.1 Integration by Parts Identities 

Integration by Parts Identities (IBP-Id's) are among the most remarkable properties 
of dimensionally regularized Feynman integrals jjj- In our case, for each of the 
integrals defined in Eq. (£[]) one can write 

where thanks to the dimensional regularization everything is well defined and the 
identity holds trivially. In the above identities the vector v M can be any of the 4 
independent vectors of the problem: ki, k 2 , pi, or p 2 , so that for each integrand 
there are 8 IBP-Id's. When evaluating explicitly the derivatives, one obtains a 
combination of integrands with a total power of the irreducible scalar products equal 
to (s — 1), s and (s + 1) and total powers fo the propagators in the denominator 
equal to (t + r) and (t + r + 1), therefore involving, besides the integrals of the family 
I t>rjS , also the families I t , r , a -u h,r+i,s and I t , r +i, s +i- 

Simplifications between reducible scalar products and propagators in the denom- 
inator may also occur, giving lower powers of the propagators. It may happen that 
some propagator disappears at all in this process; the resulting term will then give 
an integral of a simpler family (or subtopology) with (t — 1) propagators. 

As an explicit example let us consider the case of the 4-denominator topology 
of Fig. (J2J). We have three irreducible scalar products, in this topology; we choose 



0. 
0. 



(3) 
(4) 
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Figure 2: A 4-denominator topology. 



(pi ■ ki), (p 2 ■ k\) and (fci • k 2 ). Eq. (jSJ), for generic values of the indices rrii, rii and 
for a generic independent vector v^, reads: 



v» (p x ■ h) ni (p 2 ■ kxY^h ■ k 2 ) n3 



d_ 

dkf \ [kj + a] m ikj] m i(p 2 -k 2 ) 2 + a] m i(pi+p 2 - h- k 2 ) 2 + a] m ±\ 

0. (5) 



Let us take for simplicity mi = • ■ • = 1714 — 1, n\ — ■ ■ -713 — and = p^. 
Performing the derivative with respect to k\ and simplifying the reducible scalar 
products with the corresponding denominator, we write Eq. (jSJ) as follows: 








(pi ■ h 




(h ■ k 2 ) + 2 




(P2 ■ h 



(6) 



where a dot on a propagator line means that the propagator is squared and irre- 
ducible scalar products left are explicitely written. 



2.0.2 Lorentz invariance identities 

Another class of identities can be derived from the fact that the integrals I(pi), 
Eq. (PJ), are Lorentz scalars (or rather D-dimensional rotational invariant) 8\. If 
we consider an infinitesimal Lorentz transformation on the external momenta, Pi — > 
Pi + dpi, where Sp^ = e^p", and e£ is a completely antisymmetric tensor, we have 

I{Pi + 6 Pi ) = I{pi) , (7) 

Because of the antisymmetry of and because 



E 



dpn 

dp^ n 



P 
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we can write the following relation: 



Pn 



vl— 



dpi rn dp» n 



I(pi) = 



(9) 



Eq. (JUJ) can be contracted with all possible antisymmetric combination of the 
external momenta pfpj , to obtain other identities for the considered integrals. 

In our case we have two external independent momenta and we can thus con- 
struct, besides Eqs. ()3I4|) . the further identity 



(Pi • P2) Pi 



d_ 



2 u 

P2 Pi 



d 



2 u 
Pi P 2 



_d_ 
dp~i. 



/(Pi) = 0. (10) 



Let us note that, in order to obtain non-trivial identities from Eq. fjlUjl . the 
derivative with respect to the momentum pi has to be intended as a differentiation 
under the loop integral in the definition, Eq. (^), i.e. directly on the integrand of 
It,r,s- Again, this is allowed by the dimensional regularization. 

As an explicit example let us consider the same topology as in the previous 
section. If we take n\ = n<i = n 3 = and m\ = m 2 = m 3 = m 4 = 1, we have 



I(PuP2)= l{d D h}{d 1J k 2 } 



[kl + a}kl[(p2-k 2 ) 2 + a}[(p 1 +p 2 -k 1 -k 2 ) 2 + a\ 
and Eq. (fTU|) reads as follows: 








(Pi ■ h) 



(h -k 2 )-2 




+2s 



-2.1-0-- .-Q. 



2.0.3 Symmetry relations 

In general further identities among Feynman graph integrals can arise when the 
Feynman graph has some symmetry. In such a case there can be a trasformation of 
the loop momenta which does not change the value of the integral, but transforms 
the integrand in a combination of different integrands. By imposing the identity of 
the initial integral to the combination of integrals resulting from the change of loop 
momenta one obtains further identities relating integrals corresponding to a same 
graph. 
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As an example, let us consider again the topology of Fig. (|2|L with generic indices 
on the numerator and on the denominator: 



{d%}{d D k 2 } 



(pi-k 1 ) r *(jp 2 -k 1 ) n *(k 1 'k 2 ) na 



- (13) 



[kj + a\ m ^ [k%} m *l{p2-k 2 ) 2 +a] m » [( pi+p2 -k 1 -k 2 ) 2 + a\ 

The two propagators with momentum k 2 and (pi + p 2 — k\ — k 2 ) have the same 
mass. The following redefinition of the integration momentum 



h = pi + p 2 — k x — k 2 



(14) 



that consists in the interchange of the two propagators in the closed electron loop, 
does not affect of course the value of the integral; nevertheless, the explicit form of 
the integrand can change, generating non-trivial identities. 

Taking for instance n\ = n 2 = n 3 = and ni\ = m 2 = = 1, m 3 = 2 in 
Eq. f|13|) . the substitution (|14|) gives, for example, the following very simple relation: 



(15) 



Taking n x = n 2 = 0, n 3 = 1 and mi = m 2 = m 4 = 1, m 3 = 2, we get the more 
complicated idetity 






{ki ■ k 2 ) + 




(pi • fci) + 




(pa • h) 



+- 



(16) 



Summarizing, for each of the N[I t ^,s\ integrals of the family It, r ,s, Eq. (j2J), we 
have the 9 identities, Eqs. (|3I4|) and Eq. (fTUj) . involving integrals of the families up to 
It t r+i,s+i- For ( r = 0, s = 0) the number of all the integrals involved in the identities 
(for t = 6 they are 14) exceeds the number of the equations obtained (which in this 
case is 9), but when writing systematically all the equations for increasing values 
of r and s, r — 0, 1, .. , s — 0,1,..., the number of the equations grows faster 
than the number of the integrals [9 , so that at some point one deals with more 
equations than involved integrals - generating an apparently overconstrained set of 
linear equations for the unknown integrals. At this stage one can use the symmetry 
relations, somewhat reducing the number of the unknown integrals, after which one 
is left with the problem of solving the linear system of the identities. The problem 
is in principle trivial, but algebraically very lengthy, so that some organization is 
required for obtaining the solution. 
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To that aim, one can order the integrals in some lexicographic order (which means 
giving a "weight" to each integral; the weight can be almost any increasing function 
of the indices rrii, rij, such that integrals with higher indices have bigger weights) and 
then solve the system by the Gauss substitution rule by considering one by one, in 
some order, the equations of the system and using each equation for expressing the 
integral with highest weight present in that equation in terms of the other integrals 
of lower weight, and then substituting in the rest of the system. The algorithm 
is straightforward, but its execution requires of course a great amount of algebra; 
indeed, it was implemented as a chain of programs, written in the computer language 
C, which automatically runs programs written for the algebraic computer languages 
FORM ^0] and Maple [TT], reads the outputs and generates new input programs till 
all the equations are solved (and the solutions are written as a modulus of FORM 
code). 

One finds that several equations are identically satisfied (the system is only 
apparently overconstrained) , and all the appearing unknown integrals are expressed 
in terms of very few independent integrals, the Master Integrals (Mi's) for that 
family of integrals. In so doing, the resulting Mi's correspond to the integrals of 
lowest weight; but as the choice of the weight is to a large extent arbitrary, there 
is also some freedom in the choice of the integrals to pick up as actual Mi's (not in 
their number, of course!). Concerning in particular the calculation described in this 
paper, there are several cases in which two Mi's are found for a given topology or 
subtopology, while sometimes only one MI is present. It may also happen that no MI 
for the considered topology is left - i.e. all the integrals corresponding to the given 
t-propagator (sub)topology can be expressed in terms of Mi's of its subtopolgies 
with (t — 1) propagators. 

As a last remark, strictly speaking we are not able to prove that the Mi's we 
find are really the minimal set of Mi's, i.e. that they are all independent from each 
other (in the sense of the combination with polynomial factors described above); 
but in any case the number of the Mi's which we find is small, so that reducing the 
several hundred of integrals occurring in the calculation of the vertex graphs form 
factors to a few (in fact 16 Mi's, see section I27TJ) is in any case a great progress. The 
(unlikely!) discovery that one of our Mi's can be expressed as combination of the 
others would just simplify even further the calculation - without spooling, however, 
the correctness of the already obtained results. 

Concerning the number of the subtopologies, a topology with t propagators has 
(t — 1) subtopology with (t — 1) propagators, (t — l)(t — 2) subtopologies with 
(t — 2) propagators etc. It turns out, however, that most subtopologies are in fact 
equal due to symmetry relations, and the subtopologies coming from different graphs 
overlap to a great extent. For those reasons, the actual number of all the different 
subtopologies is relatively small. 

We show in Figs. dSHHJ) all the different topologies (and subtopologies; we will 
refer to them as topologies as well). In all the figures, a straight external line stands 
for an electron on the mass-shell, while the wavy external line carries the momentum 

Q. 

There are 3 independent topologies involving 6 denominators, those of Fig. Q, 
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(a) (b) (c) 

Figure 3: The set of 3 independent 6-denominator topologies present in the graphs 
of Fig. p. 



corresponding to the original graphs (a,b,c) of Fig. (JTfl. One then finds the 8 
independent topologies involving 5 denominators shown in Fig. (0J); note that the 
topologies (b,e) of Fig. (JIJ) correspond to the 6-propagator vertex graphs (d,e) of 
Fig. (JTJ), the topology (f) of Fig. (JIJ) corresponds to a vacuum polarization graph, 
(g) is a constant (i.e. does not depend on the momentum transfer Q 2 , but on a 
squared electron momentum on the mass-shell, p 2 = —a) and (h) factorizes into two 
1-loop topologies. 

Fig. (JSJ) contains all the 12 independent 4-denominator topologies; again, only 
the topologies (a,b,c,d) correspond to genuine 2-loop vertex topologies, while (e) 
corresponds to a vacuum polarization, (f,g,h) are constants and (i,j,k,l) factorize 
into two 1-loop topologies. 

Finally, Fig. (JUJ) contains all the 6 independent 3-denominator topologies; (a) 
corresponds to a vacuum polarization, the others are constants or factorizable. 

At the 2-denominator level, the only topology giving a non-vanishing contribu- 
tion corresponds to the product of two 1-loop tadpoles with squared mass a shown 
in Fig. © (g). 

2.1 The Mi's. 

For each independent topology we write the identities among the integrals of the 
associated family and solve them in terms of MPs according to the previous dis- 
cussion, see section |21 The resulting Mi's are shown in Fig. (JTJ). The Mi's are 
represented with a graph very much equal to the graph representing their topology, 
but occasional with some "decoration". When no decoration is present, the corre- 
sponding MI is nothing but the corresponding full scalar graph (first power of all 
the propagators, numerator equal to 1); that is the case, for instance, of (a), (c), (e) 
etc. of Fig. (JTJ). When a propagator line is decorated with a dot, it appears squared 
in the MI, as in the case of (d),(h) of Fig. (JTJ). Finally, the decoration can be a 
scalar product involving at least one loop momentum; the scalar product appears in 
the numerator of the integrand of the corresponding MI, as in (b), (e) etc. of Fig. 

As anticipated, there are topologies with two Mi's, as the topology (b) of Fig. 
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(a) (b) (c) (d) 




(e) (f) (g) (h) 

Figure 4: The set of 8 independent 5-denominator topologies contained in the 
graphs of Fig. ((!]). External fermion lines are put on the mass-shell p\ = p\ = —a, 
while external wavy lines carry an off-shell momentum Q = pi + P2- The topology 
(g) is evaluated on the mass-shell. 
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(a) (b) (c) (d) (e) (f) 



(g) 

Figure 6: The set of 6 independent 3-denominator (sub)topologies, (a)-(f), coming 
from the 4-denominator topologies of Fig.fjSJ) and the only non- vanishing topology 
at 2-denominator level, (g), product of two tadpoles of squared mass a. 



(JHJ) which has the two Mi's (a), (b) of Fig. 0, topologies with a single MI, as 
(a) of Fig. © which has the MI (i) of Fig. (JZj), and topologies without Mi's, i.e. 
topologies whose associated integrals can all be expressed in terms of the Mi's of 
their subtopologies; that is the case, for instance, of the topology (a) of Fig. ©. 

As a last remark, let us recall that there is some arbitrariness on the actual scalar 
integrals to be choosen as Mi's; in the case of Fig. (|7J) some of the Mi's correspond 
to graphs decorated with dots, other to graphs decorated by scalar products. The 
choice, by no means mandatory but rather somewhat accidental, was suggested by 
the convenience of later use. 

3 Calculation of the Mis. The system of differen- 
tial equations 

Once all the Mis of a given topology are obtained, the problem of their calculation 
arises. We will address the problem by the differential equations method, which 
turns out to be a really very powerful tool. The use of differential equations in one 
of the internal masses was first proposed out in ^2] > then extended to more general 
differential equations in any of Mandelstam variables in and successively used 
in ^3] for the Mi's of the sunrise diagram with arbitrary internal masses. An 
application to the 4-point functions with massless internal propagators was worked 
out in jHj for the 2-loop case and it brought to the complete evaluation of the master 
integrals for the planar [HJ H3 EH and non planar topologies [T7j . In this paper, 
we will write (and solve) the differential equations in the momentum transfer in the 
case of 3-point functions with massive fermionic propagators in QED, keeping the 
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(a) (b) (c) (d) 




(i) 0) (k) (1) 




Figure 7: The set of 16 Master Integrals (Mis). As explained in section |2~T| the 
diagrams shown are a graphical representation of the corresponding D-regularized 
integral. A dot on a propagator line means that the corresponding propagator is 
squared and an explicitly written scalar product means that the corresponding D- 
regularized integral has that scalar product in the numerator of the integrand. 
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external electron legs on the mass shell. 

Let us summarize briefly the idea of the method. To begin with, consider any 
scalar integral F(si), (we will be interested here in the Mi's, but what follows applies 
to any scalar integral as well) defined as 



F{ Si 



D i 



Ol ' ' ' Oq 



(17) 



F(si) depends in general on the three external kinematical invariants si = — p\, S2 = 
—p\ and S3 = — Q 2 = — (pi +P2) 2 , where p 2 , p\ will be later constrained on the mass 
shell p\ = p\ = —a. 

Let us construct the following quantities: 



Ojk{si) 



As F(si) depends on the Mandelstam invariants Sj, by the chain differentiation 
rule we have 



k us i 



. . d . 
Jk(si)—F(si 



(19) 



where the functions a^,jk(si) are linear combinations of the Mandelstam invariants 

Si- 

As j and k take the two values 1, 2 we obtain in that way a system of 4 linear 
equations (not all linear independent), which we can solve for the three derivatives 

d 



jj^F(si); we have in particular 







dQ- 











dp^ 



F( Sl 



where 



.4 
B 



' 1 1 
~Qi + Q 2 + 4a 
' 1 1 



Q 2 Q 2 + 4a 



(20) 

(21) 
(22) 



Assume now that F(si) is a master integral for some given topology. We can 
now substitute the right-hand side of Eq. (JT7j) in the right-hand side of Eq. (f2*Uj) 
and perform the direct differentiation of the integrand. It is clear that we obtain 
a combination of several integrals, all belonging to the same topology as -F(sj); 
therefore, we can use the solutions of the IBP and other identities for that topology 
and express everything in the r.h.s. of Eq. (|20|) in terms of the Mi's for the considered 
topology and its subtopologies. If there are several different Mi's for that topology, 
the procedure can be repeated for all the other Mi's as well. In so doing one obtains 
a system of linear differential equations in Q 2 for F(si) and the other Mi's (if any), 
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expressing their Q 2 -derivatives in terms of the Mi's of the considered topology and of 
its subtopologies; due to the presence of the Mi's of the subtopologies the equations 
are in general non-homogeneous. 

At this point we can impose the mass-shell conditions, and the general structure 
of the system reads 

^M,(Aa,g 2 ) = Aj(D, a, Q 2 )Mj(D, a, Q 2 ) 

j 

+ B h {D, a, Q 2 )N k (D, a, Q 2 ) , (23) 

k 

where the Mi(D,a,Q 2 ) are the Mi's of the topology with the electron legs on the 
mass-shell, Nk(D, a,Q 2 ) the Mi's of the subtopologies, we have made explicit the 
dipendence on D for later use (note that the above equations are exact in D) and the 
coefficients Aj(D,a,Q 2 ), Bk(D,a,Q 2 ) are rational factors depending on D,Q 2 and 
the electron squared mass a. As will be apparent by the examples, the singularities 
of Aj(D, a, Q 2 ), B k (D, a, Q 2 ) in the variable Q 2 , such as 1/{Q 2 + 4a) and l/Q 2 , and 
correspond to the thresholds and pseudothresholds of the corresponding Feynman 
graphs. 

It is clear that the procedure can be repeated in principle for the other Man- 
delstam variables as well. We are not interested in this further equations, as in the 
case we are considering the other Mandelstam variables are the invariant masses of 
the electrons, which we keep frozen on the mass-shell. 

As already observed, the Eqs. (|2^j) for the Mi's Mi(D, a, Q 2 ) of a given topology 
are not homogeneous, as they may involve the Mi's Nk(D, a, Q 2 ) of the subtopolo- 
gies. It is therefore natural to proceed bottom up, starting from the equations for 
the Mi's of the simplest topologies (i.e. with less denominators), solving those equa- 
tions and using the results within the equations for the Mi's of the more complicated 
topologies with additional propagators, whose non-homogeneous part can then be 
considered as known. 

3.0.1 The boundary conditions 

Some comments on boundary conditions. As already observed, the coefficients of 
the differential equations Eq. (J23j) are in general singular at Q 2 = and Q 2 = —4a; 
correspondingly, the solutions of the equations can have a singular behaviour in those 
points. But we know that the Vertex integrals are regular in Q 2 at Q 2 = when the 
electron lines are on the mass shell, a qualitative result which can easily verified, 
when needed, by direct inspection of the very definition of the amplitudes as loop 
integrals. It turns out that the qualitative information provided by the regularity 
behaviour (implying the absence, in the Q 2 —* limit, of terms in l/Q 2 or \nQ 2 ) is 
completely sufficient for the quantitative determination of the otherwise arbitrary 
integration constants which naturally arise when solving a system of differential 
equations. 
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3.1 Laurent series expansion in e 

The system of differential equations Eq. is exact in D, but we are interested, 
in any case, on the Laurent expansion of the solutions in powers of e = (4 — D)/2. 
It turns out that our 2-loop integrals have up to a double pole in e (which can be of 
ultraviolet or infrared origin), so that quite in general we will expand the two loop 
Mi's as 

n 

Mi(D,a, Q 2 ) = £ e>M«\a,Q 2 ) + 0(e^) , (24) 
j=-2 

where n is the required order in e. We will present for all the considered integrals 
the coefficients of the e expansion up to the zeroth order included. In some cases, 
however, we had to expand intermediate results up to the term of fourth order in 
e. That depends on the fact that some of the Mi's, which appears in the non- 
homogeneous part of a system of differential equations for more complicated Mi's, 
can be multiplied by coefficients which are also singular in e. 

When expanding systematically in e all the Mi's (including those appearing 
in the non-homogeneous part) and all the D-dependent coefficients of Eq. ()23|). 
one obtains a system of chained equations formed by the equations valid for each 
power of e. The first equation corresponds to the coefficient of double pole in e 
of the equation, and involves only the coefficients M.\ 2 \a,Q 2 ) as unknown; the 
next equation, corresponding to the single pole in e, involves the M-~ 1 \a,Q 2 ) as 
unknown, but can in general involve M^~ 2 \a, Q 2 ) if some of the coefficients contains 
a power of e; such a term in M.\ 2 \a,Q 2 ) can be considered as known once the 
equation for the double pole has been solved. For the subsequent equations we have 
the same structure: they involves the coefficient M^'(a,Q 2 ) at the order j in e as 
unknown and the previous coefficients as known non- homogeneous terms. 

Let us note that the homogeneous part of all the equations arising from the e 
expansion of Eq. ()23j) is always the same and obviously identical to the homogeneous 
part of Eq. (}2"Hj) at D = 4, i.e. e = 0. It is natural to look for the solutions of the 
chained equations by means of the Euler's method of the variation of the constants, 
using repeatedly the solutions of the homogeneous equation, as we will show in some 
examples in the next section. 

General algorithms for the solution of the homogeneous equations are not avail- 
able; it turns out however that in all the cases considered in this paper the homoge- 
neous equations at D = 4 have almost trivial solutions, so that Euler's formula can 
immediately be written. With the change of variable 

x = —, v —= , (25) 

y/Q 2 +4a+y/Q* 

all integrations can further be carried out in closed analytic form, the result being 
a combination of the Harmonic Polylogarithms introduced in |Hj (see also [6 j for 
their numerical evaluation), a generalization of the already widely used Nielsen's 
Polylogarithms [21 E]- 
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As a last remark some comments on the arbitrariness of the choice of the Mi's. 
For the topologies we considered, we had at most 2 Mi's, which means that we had 
to solve in principle a linear system of two equations or an equivalent second-order 
differential equation. However, the freedom in the explicit choice of the Mi's can 
play an essential role in simplifying the calculation. It turns out, in fact, that if we 
choose the two Mi's with a different leading singular behaviour in e, the system of 
the two coupled first-order linear differential equations does in fact decouple. 5 As 
a result, instead of solving a second-order differential equation we can solve simply 
two first-order equations. We will show how to exploit this possibility in the solution 
of the systems considered in sections 14.21 and 14.31 

4 Explicit calculations 

In this section the equations for three topologies and their solutions are discussed 
in some details. We chose a 4-, a 5- and a 6-dominator topology, shown respectively 
in Fig. © (a), Fig. (JHJ) (a) and Fig. (jHJ) (b), to illustrate the algorithm for the 
solution of the corresponding system of differential equations. 

The 4-denominator topology is the simpler among the three cases, since it has 
only one MI. Correspondingly, the system in Eq. ()23|) reduces to a single first-order 
linear differential equation - whose solution is therefore trivial. 

The other two topologies are more difficult. They have both two Mi's and 
therefore, in both cases, we must solve in principle a system of two first-order coupled 
linear differential equations. As we have already remarked in section |3~T1 the explicit 
form of the system depends on the choice of the pair of Mi's. Indeed, we choose 
in both cases two Mi's with different leading behaviour in e, such that the system 
decouples order by order in e. 

Before to go on, two remarks have to be done. 

The first one is the following. We are interested in really 2-loop diagrams, but, as 
we have seen, in the reduction to the Mi's we encountered topologies which factorize 
in the product of two 1-loop topologies. In particular, two Mi's have this structure; 
they are shown in Fig. (J7|) (k) and (p) and they consist respectively in a product of 
two bubbles in the Q-channel and of one bubble in the Q-channel and a Tadpole. 
The same algorithm explained in sections [2HH1 wa s applied, therefore, to the 1-loop 
problem and the results are shown in appendix [Bl where we discuss in particular 
the solution of the differential equation for the bubble with two massive propagators 
with equal squared mass a. 

The second remark concerns the 3-denominator Mi's. They constitute the sim- 
plest non-trivial 2-loop topologies of the entire pyramid of Mi's, and thus they are 
present in the non-homogeneous part of the systems of differential equations for 
all the other topologies. For the two Mi's of the Sunrise with two equal-mass and 

5 The decoupling can be, in some cases, exact in e, which means that a combination of integrals 
diagonalize the system without expanding it in powers of e. More in general, the homogeneous 
equation at D = 4 is brought to acquire a triangular form, with the first differential equation which 
contains only one of the Mi's and the second equation which involves both Mi's. 
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one mass-less propagator, Fig. (J7J) (m) and (n), a system of two non- homogeneous 
first-order differential equations can be established, the non-homogeneous part con- 
sisting essentially on the product of two massive 1-loop Tadpoles, times a ratio of 
polynoms in Q 2 , a and e. But, if we consider the Sunrise with two mass- less prop- 
agators, analogous to that one in Fig. (J7J) (o), but with the external leg off-shell, 
this is no longer possible. The resulting system is, in fact, homogeneous, as we 
can understand from the fact that contracting a propagator line we have at least a 
product with a mass-less Tadpole, which vanishes in dimensional regularization. In 
this case the conditions of regularity of the integrals in Q 2 = are not sufficient to 
determine the boundary conditions, as explained in section 13.0.11 In this situation 
we are forced to evaluate the integrals by direct integration, as we did for the MI in 
Fig. (J7J) (o); but that is not a problem, given the simplicity of the integrals. 



4.1 The full calculation for the topology in Fig. (03) (a) 

The topology in Fig. (j5J) (a) has only one MI. We choose the simpler one, i.e. the 
scalar integral itself, Fig. (J7|) (e): 



F{e,a,Q 2 




, 2(4-D) 



{d D k 1 }{d u k 2 } 



1 



(26) 



where the exlicit expressions of the denominators T>i are given in appendix [X] 

The first-order linear differential equation in the variable Q 2 , which we obtain 
by the methods described in the previous sections, reads 



dF(e,a,Q* 
dQ 2 



1 

g - 3e) 
4a 



(l-4e) 
(Q 2 + 4a) 

1 

~Q 2 ~ ~ 



F(e,a,Q 2 



(Q 2 + 4a) 



(27) 



where the 3-denominator diagram on the non-homogeneous part of the equation is 
the MI of Fig. ((7j) (o), function only of the squared mass a; its expansion in Laurent 
series of e is given in Eqs. (J85H87J) . 

As we can notice, from Eq. (}2Tj) we see that its solutions can be singular at 
Q 2 = —4a and Q 2 = 0. The integral we are considering, Eq. (J2T?|) . is indeed singular 
at the physical threshold Q 2 = —4a, but regular at the pseudothreshold Q 2 = 0. 
This allows us to get the initial condition directely from the differential equation. 
In fact, multiplying Eq. (|27|) for Q 2 and taking the limit Q 2 — > 0, the left-hand side 
simply vanishes; the right-hand side gives us: 



F(e, a, Q 2 = 0) = -0- (28) 

Even if in this particular case it is possible to find a solution of Eq. (|27J) exact 
in D = 4 — 2e, jTH], in this section we look for a solution expanded in Laurent series 
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of e: 



o 



F(e, a, Q 2 ) = ^e^a, Q 2 )+ 0(e) 



(29) 



i=-2 

The homogeneous equation at D = 4, i.e. e = 0, is 



df(a,y) 
dy 



2 Li/ (|/ + 4a) 



1 1 1 



/(o,y) > 



(30) 



whose solution is 



f(a,y) = kJl + 



4a 



(31) 



y 



where k is a normalization constant. 

Order by order in e, we obtain the solution of the non- homogeneous equation by 
means of the method of the variation of the constant k (Euler's method). Substi- 
tuting, into Eq. (J27J), Eq. (j2HJ) and the expansion in e of the 3-denominator integral 
on the non- homogeneous part, Eq. (jH^j) . the result reads as follows: 



where the explicit values of the constants Ci are given in Eqs. (JHSHSZI) ■ 

The determination of the constants ki is made by imposing that the solution, Eq. 

(I32j) satisfies the initial condition, Eq. (J28|) . or, which is the same thing, imposing 

the regularity of the solution at Q 2 = 0. 

It is then useful to express the result in terms of the variable x, defined in Eq. 

(I25J) . The explicit form of the solution up to the zeroth order in e is given in Eqs. 

(|105H107J)- 

4.2 The full calculation for the topology in Fig. (JH) (a) 

The topology in Fig. Q (a) has two Mi's. In order to decouple the system of 
differential equations order by order in e, we choose the couple of Mi's given by Fig. 
(J7J) (c) and (d), i.e. the fully scalar integral and the scalar integral with a squared 
propagator: 




(32) 




(33) 



(34) 
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The corresponding system of first-order linear differential equations in the vari- 
able Q 2 reads 

dF^a.Q 2 ) 1 



dQ 2 



dF 2 (e,a,Q 2 ) 
dQ 2 



1 (l~2e) 
g 2 + (Q 2 + 4a) 



Fi(e,a,Q 2 ) 



F 2 (e ) a,g 2 ) + Q«(e,a,g 2 



e 

2a 



Q 2 (Q 2 + 4a) 



F!(e,a,Q 2 ) 



1 + (l + 2e) 



F 2 (e,a,g 2 ) + fi( 2 )(e,a,g 2 ), 



g 2 (g 2 + 4a)_ 

where the functions fiW(e, a, g 2 ) are the following combinations of Mi's: 

I f (l-2e) 2 (3-4e) 



(35) 



(36) 



n (1 \e,a,Q 2 ) 



1 (l-3e)(l-2e) 
2ag 2 (I - 4e) 

(3 - 20e 2 + I6e 3 ) 
~16a 2 (l-4e)(g 2 + 4a) 




2e [ I6a 2 (l -4e)g 2 

(II-I08e + 256e 2 -I68e 3 ) 
4a(l -4e)(g 2 + 4a) 2 



3(5- I6e+ I2e 2 



(g 2 + 4a) 3 
3(l-2e) 2 (l-e) 



J ^Xy^ e \ l 6a 3(! _ 4e )Q2 



3(l-e-4e 



4f 



-3^ 



16a 3 (1 - 4e)(g 2 + 4a) 4a 2 (1 - 4e)(g 2 + 4a) 2 
9(1 -3e + 2e 2 



a(g 2 + 4a) 3 
(3 - 5e + 2e 2 




(k k) f ( 3 ~ 5e + 2£2 ) 
^ ^ "\l6a3(l-46)g 2 

(2 - 7e + 5e 2 ) 



16a 3 (l - 4e)(g 2 + 4a) 4a 2 (l - 4e)(g 2 + 4a)' 

(3-9e + 6e 2 )l 2 

T (e,a) , 



fi (2) (e,a,g 2 ) 



ae(g 2 + 4a) 3 

e(l-2e)(l-3e) 
2a 2 (1 -4e) 

(l-2e)(l-3e) 



(37) 



4a 



g 2 (g 2 + 4a) 
1 1 




1 



(g 2 +4a) 2 



1 




4ag 2 4a(g 2 + 4a) 

(l-2e)(l-3e)(2-3e) 



+ 



1 



4a(g 2 + 4a) (Q 2 + 4af 



(l-e)(l-2e) 



8a 2 (l -4e) 



-e- 



4ag 2 



4a 2 



4ag 2 4a(g 2 ^a) (g 2 +4a) : 



O 2 
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(l-2e)(l-lle+12e 2 



32a 3 (l-4e) 



1 

~Q 2 



(1 - 26) 
4e 



(3-4e)(l-4e) 
16a 2 Q 2 



(Q 2 + 4a)J 
(3 - 36e + 296e 2 



288e 3 



16a 2 (l -4e)(Q 2 + 4a) 2 



3(1 - 16e + 20e 2 ) 3(5 - 6e)(l - 6e)l \ r J^\ 
"2a(l -4e)(Q 2 + 4a) 3 + (Q 2 + 4a) 4 J jN^T' 



3(l-e)(l-2e)(l + 4e) 



64a 4 (l-4e) 



1 

Q 2 



l 



3(l-e)(l-2e) 
2ae 

(l-7e + 36e 2 ) 
'2a(l-4e)(Q 2 + 4a) 3 

(l-e)(3-31e + 122e 2 



(Q 2 
(1 



4a) _ 

6e + 24e 2 ) 



(1 - 46) 

16a 2 g 2 16a 2 (l -4e)(Q 2 + 4a) 

3(1-66) ]\ s n. 



+ 



(Q 2 + 4a) 4 
- 104e 3 



(1 



2ae 



128a 4 (l-4e) 

e(l~4e) , (3 
16a 2 g 2 



1 

32e + 130e 2 



1 



(h ■ h) 



(Q 2 + 4a) 
120e 3 ) 



16a 2 (l -4e)(Q 2 + 4a)^ 



(3 - 29e+ 134e 2 -168e 3 ) 3(1 - 2e) (1 - 6e) 



2a(l-4e)(Q 2 +4a) J 



(Q 2 + 4a) 4 



T 



e, a 



(38) 



where T(e, a) stands for the tadpole explicitly definited in Eq. (fTUj) . 

For what initial conditions are concerned, we know that F\(e, a, Q 2 ) is analytic 
in Q 2 = 0. For Euclidean momenta, the limit Q 2 — > (which implies Q — > 0) can 
be recovered by the limit P2 — > — Pi- Taking this limit directly within the integrand 
of Eq. (jnHJ) we obtain 



F 1 (e,a,Q 2 = 0) 



(2-3e)(l-3e) 
2a 2 (l -4e) 



(l-6) 2 (l-2e) 
2a 3 (l - 2e) 



T 2 (e,a). (39) 



Once we have the initial condition for -F\(e, a, Q 2 ), we can calculate the initial 
condition for the second MI, F 2 (e, a, Q 2 ), directly from Eq. (|33j) . In fact, because of 
the analyticity of Fi(e, a,Q 2 ) in Q 2 = 0, we can multiply Eq. (|33j) by Q 2 and take 
the limit Q 2 — > 0. The left-hand side vanishes and we find the following relation: 



F 2 (e,a,Q 2 = 0) 



3e(l-3e)(2-3e) 



8a 3 (1 -4e) 

(16e 4 - 16e 3 - 76e 2 + 226e - 171) 
+ 32a 4 



T 2 (e,a). (40) 
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We look for a solution of the system of Eqs. 
of the Laurent series in e: 



in terms of the coefficients 



F x (e, a, Q 2 ) 
F 2 (e,a,Q 2 ) 



e l F^(a,Q 2 ) 



i=-2 



O(e) 



O(e) 



(41) 
(42) 

, the systematic ex- 



As immediately seen by direct inspection of Eqs. 
pansion in powers of e gives a triangular system in which the second equation, at 
order % in the e expansion, consists of a homogeneous part with F^ 2 \a ) Q 2 ) only, 
and a non-homogeneous part which contains expansion terms of order lower than 
i. The first equation for F± (a,Q 2 ), on the contrary, contains both F- 1 \a,Q 2 ) and 
F^ 2 \a,Q 2 ) in the homogeneous part; but as a first step we can solve the second 
equation for F^ 2 \a, Q 2 ), substitute the result in the first equation and split again 
the first equation in a new homogeneous part, which now involves only (a, Q 2 ), 
and a non- homogeneous part which is known. As a result, the original system splits 
into two homogeneous decoupled equations, which are 



dfi(a,y) 
dy 

dy 



111 1 

- — h ~, r fi(a>,v) 

2 [y (y + 4a)\ Jn ,y) 

- + 7 } A \ /a(a,J/), 
1/ (2/ + 4a)_ 



with solutions 



h(a,y) 

f2(a,y) 



h 

^y(y + 4a) ' 
fc 2 

y(y + 4a) ' 



(43) 
(44) 

(45) 
(46) 



By means of the Euler method we can find, order by order in e, the solution 
of the non-homogeneous system, solving the two first order differential equations in 
the order of the two quadrature formulas: 



F?\a,Q 2 ) 



Q 2 (Q 2 + 4a) 



Q 



dyy{y + 4a) 



2a \y (y + 4a) 



F?\a,Q 2 ) 



(y + 4a) 
1 



,(2) 



(47) 



4a 



dyVy(y + ^ a ) 



i 



(i) 



d/ + 4a) F ^ 1 



a,y) 



(4£ 
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where, for simplicity, we put: 

o 

Q^(e,a,Q 2 ) = ^6^ (1) (a,g 2 ) + 0(e) 

i=-2 


rt 2 \e,a,Q 2 ) = ^e^ (2) (a,Q 2 ) + 0(e) 



(49) 
(50) 



i=-2 



The determination of the constants k[ and kf* 1 is made imposing the initial 
conditions Eqs. (pM l4*U|) . The solution, expressed in terms of the variable x, is given 
in Eqs. (113711131)1) . 

4.3 The full calculation for the topology in Fig. (03) (b) 

The topology in Fig. (jSJ) (b) has two Mi's. We choose the Mi's corresponding to 
Fig. (J7J) (a) and (b), i.e. the fully scalar integral and the scalar integral with the 
scalar product (ki ■ k 2 ) in the numerator of the integrand. Also in this case the 
choice of the Mi's diagonalizes the system in the limit e — > 0. 



F 2 (e,a,Q 2 




/'o 



2{4 - D) l{d%}{d%} 



1 



v 1 v 2 v 9 v u v 14 v 



15 




(h-k 2 )=^~ D) {d D h}{d D k 2 } 



h ■ k 2 



V 1 V 2 V 9 V 11 V U V 15 



The system of first-order linear differential equations is the following: 



(51) 



(52) 



dF^e, a, Q 2 
dQ 2 



dF 2 (e,a } Q 2 ) 
dQ 2 



;i + 26) 



1 (l-2e) 



2e 
a 



[Q 2 (Q 2 + 4a)J 
(l-2e1 



F x (e, a, Q 2 



[Q 2 (Q 2 + 4a)\ 
eF^a^Q 2 ) 
(I - 2e) 



F 2 {e,a,Q 2 ) + Q^{e,a,Q 2 



[Q 2 (Q 2 +4a)J 



where the functions f2^(e,a, Q 2 ) are defined as follows: 



fi (1) (e,a,Q 2 



e 
a 



Q 2 (Q 2 + 4a) 




3(1 -2e) 



4a 1 
Q^'O 2 



(Q 2 + 4a) 




3(2 - 3e) 
4a 3 



4a 1 
Q^'Q 2 



(53) 



F 2 (e,a,Q 2 )+rt 2 \e,a,Q 2 ), (54) 
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(Q 2 + 4a) 




5 

+ 2a 



1 



3(l-2e)(l-3e) 
4a 2 (1 -4e) 



(l-2e)(l-3e) 



Q 2 (Q 2 + 4a) 
1 




Q 2 (g 2 + 4a) 




4a 



16a 2 e [Q 2 (Q 2 + 4a) (Q 2 + 4a) 2 



(17- 116e + 248e 2 - 128e 3 ) 




Q 2 (Q 2 + 4a) 



64a 3 e(l - 4e) 
(3-4e)(l-4e)(l-22e) 1 
64a 2 e 2 g 1 
(3 - 162e + 968e 2 - 2240e 3 + 1536e 4 ) 1 

64a 2 e 2 (l - 4e) (Q 2 + 4a) 2 

3(1 - 2e)(l + 34e- 140e 2 + 120e 3 1 

8ae 2 (l -4e) (Q 2 + 4a) 3 

3(l-2e)(5-6e)(l-6e) 1 1 /fX 
4^ (Q 2 + 4a) 4 J 



+ 



Q 2 (Q 2 + 4a) 



(Q 2 + 4a) 4 J 

3(l-e)(9-2e-112e 2 
64a 4 e(l - 4e) 
3(1 -e)(l -4e)(l - 22e) 1 
32a 3 e 2 g 4 
3(1 -e)(l -48e + 196e 2 - 128e 3 ) 1 

32a 3 e 2 (l -4e) (Q 2 + 4a) 2 

3(1 -e)(l -2e)(l - 37e + 96e 2 ) 1 

4a 2 e 2 (l-4e) (g 2 + 4a) 3 

9(l-e)(l-2e)(l-6e) 1 1 /^v ( ^ 



2ae 2 



(Q 2 +4a) 4 J 



fc 2 ) 



(l-2e)(l-3e)(2-3e) 



32a 3 e(l - 4e) 



g 2 (g 2 + 4a) 



4a 



(g 2 +4a) s 



3(2-3e) 
16a 3 



4a 1 1 

f 



g 4 g 2 (g 2 +4a). 
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X- 



3(2-3e) 
32a 3 



4a 



g 4 Q 2 (g 2 +4a). 



O 2 



(l-e)(3-89e+436e 2 -576e 3 + 16e 4 ) [ 1 
128a 4 e 2 (l-2e)(l-4e) 
(1 -e)(l -38e+ 100e 2 ) 1 



g 2 (g 2 +4a)J 



32a 3 e(l - 2e) g 4 
(l-e)(3-84e + 310e 2 - 208e 3 ) 1 

32a 3 e 2 (l - 4e) (g 2 + 4a) 2 

(l-e)(3-89e + 374e 2 - 408e 3 ) 1 



8a 2 e 2 (l-4e) 
3(l-e)(l-2e)(l-6e) 



(g 2 + 4a) ; 



2ae 2 



(g 2 + 4a) 4 



T 2 (e,a), 




+ 



(l-2e)(l-4e) 



2e 



4a 



+ 



g 4 (g 2 + 4a) 



(1 - 2e) 1 

8ae g 2 



(1 -2e)(2-3e) 
4a 2 e 



(1 - 3e) 




+ 



(2-3e)(l-4e) 1 

m g 4 



g 2 (g 2 + 4a) 




(P2 • fcl) 



a (g 2 + 4a) 




(1-26) 1 
2e g 2 



(1-66) 



2e (g 2 +4a) 



V (l-2e)(l-3e) 1 

2a(l-4e) g 2 




+ 
+ 
+ 



(l-2e)(3-4e)(l-4e) 1 

1606 2 g 4 

(3 - 18e + 4e 2 + 136e 3 - 160e 4 ) 1 
64a 2 e 2 (l - 4e) g 2 
3(l-2e)(5-6e)(l-6e) 1 

4^ (g 2 + 4a) 3 

(3 - 74e + 410e 2 - 816e 3 + 504e 4 ) 



8ae 2 (l -4e) 



(g 2 + 4a)^ 



(3 - 18e + 36e 2 + 40e 3 - 96e 4 ) 1 



64a 2 e 2 (l-4e) (Q 2 + 4a) J 

3(1 - e) | (1 - 2e)(l - 4e) 1 (1 - 32e 2 + 72e 3 ) f 1 



2ae 2 



4a 



g 4 



16a 2 (l-4e) Vg 2 
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(Q 2 + 4a) 
(1 -5e + 4e 2 + 12e 3 



3(1 -2e)(l -6e) 



1 



(Q 2 + 4a) 3 



l^A^-O-^ {ki ■ 



2a(l - 4e) (Q 2 + 4a) 2 J 



feu) 



(l-2e)(2-3e)(l-3e) 
16a 2 e(l - 4e) 



(2-3e)(l-4e) 1 (l-2e)(2-3e) 



8ae 



(Q 2 + 4a) 2 



Q 4 

e 



1 

Q~ 2 



8ae 

(l-e)(l-4e) 1 



2ae 



Q 4 



(l-€)(l-2e) 
8a 2 e 



O 2 



Q 2 (Q 2 + 4a) 2 _ 

(1-e) / (l-4e) 1 (3-25e-71e 2 + 58e 3 + 32e 4 ) 
2a \4a(l-2e)Q4 + 16a 2 e 2 (l-2e)(l-4e) 
1 \ 3(l-2e)(l-6e) 1 



1 

~Q 2 



(Q 2 + 4a), 
3(1 - 7e + 14e 2 ) 



4ae 2 (l-4e) (Q 2 + 4a) 2 



(Q 2 + 4a) J 



T 2 (e,al 



(56) 



For what concerns initial conditions, we know that Fi(e, a, Q 2 ) is analytic in 
Q 2 = 0. For Euclidean momenta, the limit Q 2 — > (which implies Q — > 0) can be 
recovered by the limit p 2 — > —pi- Taking this limit directly within the integrand of 
Eq. f)5ip we have: 



F 1 (e,a,Q 2 = 0) 



■<D- 

3e(2-3e)(l-3e) 
~4a 3 (l-4e)(l + 2e)" 




3(2-3e)(l-3e) 



64a 3 e 



(1 - e) 2 (9 + 3e - 160e 2 - 196e 3 , 2 



T z (e,a). (57) 
64a 4 e(l-2e)(l + 2e) V ' ; V ' 

We can find the initial condition for F 2 from Eq. ([53)1 . multiplying by Q 2 and 
taking the limit Q 2 — ► 0, or performing the limit directely inside the integral of Eq. 
(J52j) . What we get is the following expression: 



F 2 (e,a,Q 2 = 0) 



(2-3e)(l-3e) 

8^6 




+ 



(2-3e)(l -3e) 
32a 2 e 2 



e 
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We look for a solution of the system of Eqs. (jSnHSU) expanded in Laurent series 



of e: 



Ffaatf) = J]e i J F; (1) (a,Q 2 ) + 0(e) 



i=-2 




F 2 (e,a,Q 2 ) = £ e^ (2) (a, Q 2 ) + O (e 



(59) 



(60) 



i=-2 



According to the previous remarks, the solution can be built, order by order, by 
the method of the variation of the constants of the associated homogeneous system, 
which reads 



dfi{a,y) 
dy 

df^hy) 
dy 



1 1 

+ 



y (y + 4a) 

l l 

+ 



h{a,y) 



(61) 



(62) 



y (y + 4a)_ 

As we can see, the homogeneous system is completely diagonalized in the limit 
e — ► 0. The solution of the system is the following: 



h(a,y) 



y{y + 4a) ' 

k 2 



(63) 
(64) 



By means of the Euler's method we can find, order by order in e, the solution of 
the non-homogeneous system: 



Q 2 (Q 2 + 4a) 



F}%(a,y) + W(a,y) 



(y + 4a) 
4 

a(y + 4a) 
1 

7§W+4aT 
1 1 



dyy{y + 4:a) 
2 (\ 



, 1 1 
-2 ( - + 



y (y + 4a) 



a \y (y + 4a) 



-A; 



(1) 



(65) 



V£5/ + 4a) F£\(a,y) 



y (y + 4a) 
where, for simplicity, we put: 



fS(«,2/) + nr(M) 



(2), 



-k) 



(2) 



fi«(e,a,g 2 ) = ^e^ (1) (a,Q 2 )+0( e ) 



i=-2 




fi( 2 )(e,a,Q 2 ) = ^e^ (2) (a,Q 2 ) + 0(e) 



(66) 

(67) 
(68) 



i=-2 
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The determination of the constants fc^ and kf* 1 is made imposing the initial 
conditions Eqs. (|57 | 158 )) . The solution, expressed in terms of the variable x, is given 
in Eqs. (1144111461) . 



5 Results for the Mi's 

We list in this section all the Mi's necessary for the calculation of the 2-loop vertex 
diagrams of Fig. ((TJ) . 

We give them as a Laurent series in e = (4 — D)/2 and we express the coeffi- 
cients of the series in terms of harmonic polylogarithms of the variable x, already 
introduced in Eq. (J2HJ) 

VQ 2 + 4a - v/Q2 

x - 



For brevity we present the results only up to the zeroth order in e, but of course 
the method allows us to calculate any order in e. All the coefficients of the e expan- 
sion depend of course on a and Q 2 (or the above variable x) ; for reducing the size of 
the formulas, we will not write anymore the dependence of the coefficients on those 
variables. The analytic results are expressed in terms of Harmonic Polylogarithms 
of argument x. Definition, notation and properties of the Harmonic Polylogarithms 
can be found in jSJ |H] . 

The scale /xo is the regularization scale and a = m 2 e is the only mass involved (in 
our case the mass of the electron). 

The explicit values of the T>i is given in appendix 1X1 

For what concerns the normalization of our integrals, we define the 1-loop Tad- 
pole with mass a as 

T{e,a)=^- D) j{d D k}^-. (69) 
We further define {d D k} in order to have 6 : 



T(e,a)= - -, (70) 

so that 

rl D k 

w = D) • (71) 

with D = 4 - 2e, e = (4 - D) /2. 

All the results can be downloaded as an input file for FORM in [T8j . 



6 With this normalization the tadpole T(e, a) of the present paper is 4 times the corresponding 
quantity of [T5 ] . [H ] and 
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5.1 Topologies with t = 3 



1 



v 2 v e v w 



i=-2 



rs^^^, [ki ■ k 2 



2(4-D) I t A D 



{d n k 1 }{d D k 2 } 



D, i h- h 



T> 2 V 6 V W 



-2e 



(72) 
(73) 
(74) 
(75) 



i=-2 



As already said above, from now on we write for short Ai, Bi instead of Ai(a, Q 2 ) 
and Bi(a,Q 2 ). Referring to ^H] for more details, we find 



A. 
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A-x 
a 
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a 
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a 2 
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11 13 
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(77) 



1 
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1 
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1 
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1 1 
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1 
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X 



1 

24 


" 1 2 

x z 


11 
~ 24 


"1 

— V x 

X 


11 
12 ' 


1 

~ 48 


13 




2 )-ll 



(78) 
(79) 
(80) 



x 



X 



(1-x) 



1 

12 



- x 2 - 7 



X 



H(0,x) + 



H(0,x) 
2 



x (1 — x) 



+ 



(1 -x) 2 _ 



#(0,0, x). 



51) 




2(4-r>) 

P0 



i 



VxV 2 V 15 



32) 

S3) 
34) 
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As mentioned in section we calculated this MI directely by means of Feynman 
parameters. We found: 



a 

cu 

a 

Co 



1 

2 
5 

4 
11 



-5i = ---2C(2) 



(7 



55) 
56) 
57) 



where: 



a 
a 

a 



-1, 
-3+ I 

-7- 1 



/if - D) / 



V 6 V 7 V 12 



i=-2 



(1-x) 
2 



(1-x) J 
2#(-l,0,x 



if(0,x), 
C(2) - 3if(0, ar) - #(0, 0, x) 



59) 



(90) 
(91) 



(92) 



where 



e 



v 6 v 7 v 15 



/'o 



-2e 
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~2 ' 




17 


a 


~T 


Fo 


59 
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(93) 
(94) 



(95) 
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(97) 
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Topologies with t = 4 
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1 - 



#(0,x) 
2 
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-2e 
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where: 



J-2(X) 

J-i(x) 
J (x) 



1 

2' 
5 

2' 
19 



+ 2 



(l-x) (l + x) 



^— ]{C(2)tf(0,x) + #(0,0,0, x)} 



-H(0,0,x). 



(110) 
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(112) 
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aL^ 2 
aL_i 

aL Q 
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5.3 Topologies with t = 5 
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5.4 Topologies with t = 6 
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+4ff(0, 1, 1, 0, x) - 8H(1, 0, -1, 0, x) + 8i7(l, 0, 0, 0, x) 
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6 The 6-denominator reducible diagrams 

The other diagrams of Fig. ((TJ) are all reducible diagrams. We give in this section 
their result. 
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7 Expansion for Q 2 ^> a 

We list, in this section, the asymptotic expansion of the 6-denominator vertex dia- 
grams given in the previous sections, in order to show their behaviour for momentum 
transfer larger than the mass. 

Putting y = Q 2 /a , L = \ny and keeping terms up to the order (1/y) 5 , we have: 
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+4((2)L 2 - 4L 2 + 4L 3 - 2L 4 , 
296 

+ 18C(2) + — C 2 (2) + 15C(3) - 20 + 18C(2)L - 16C(3)L 
o 

70 

+3LL - 16C(2)L 2 + 37L 2 - —L 3 + 8L 4 , 
724 w . 1184 >9 . . 248 w . 208 w . 

-— c(2) - — c 2 (2) - — c(3) + iss - —cm 



2204 

+64C(3)L-^ r L + 64C(2)L 2 
-32L 4 , 

-ic 2 (2)-2C(3)L-l/A 



1948 



-L + 112L 3 



(167) 
(168) 



(169) 
(170) 
(171) 
(172) 
(173) 

(174) 

(175) 



176) 

;i77) 
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-2 - 2C(2) + l( 2 (2) - 4C(3) + 4C(2)L + 4C(3)L - 2L 

T + Y C(2) " ^ 2(2) + 14C(3) + T L " 14C(2)L 

-12C(3)L + ^L 2 -^L 3 -iL 4 , (179) 

_ isl 2 + —L 3 + -L 4 , (180) 
54 9 6' V ; 

146447 25325 , , >2 .„, 533 >/0 , 52955 r 

loses + — c(2) " 14C (2) + -« 3) + 16T L 

-^C(2)L - 140C(3)L + "-fl? - r §L> - f 2 L' . (181) 



-2e 5 



(183) 



L 2 
2 ' 

2L - 2L 2 , (184) 

2-11L + 8L 2 , (185) 
152 

-14 + — L-32L 2 , (186) 
o 

-C(3) + 2C(2)L-iL 3 , (187) 

4C(2) + 4C(3) -2L- 8C(2)L - 3L 2 + , (188) 

7 37 8 

-3 - 22C(2) - 16C(3) + -L + 32((2)L + -L 2 - -L 3 , (189) 

47 304 70 

- + ^-C(2) + 64C(3) + ^L- 128C(2)L - 92L 2 

+ f^ 3 , (190) 
f( 2 (2) + a3)L+ 7 -(m 2 + ^L\ (191) 
-8 - 2C(2) - yC 2 (2) + 6C(3) - 2L + 6((2)L - 4((3)L 

-14C(2)L 2 + L 2 + 7 -L 3 - h\ (192) 

93 1Q9 8Q 

48 + yC(2) + ^-C 2 (2) - 37C(3) + ^L- 25((2)L + 16C(3)L 
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+ ^L 2 + 56C(2)L 2 -^L 3 + ^L 4 , (193) 

36349 362 w . 768 >2/ . J>/n . 3395 r 

-C(2) - — C 2 (2) + 184C(3) - —L + 88C(2)L 



162 9 " v 7 5 ' v ' " v 7 27 
-64C(3)L - ^L 2 - 224C(2)L 2 + ^L 3 - f L * . (194) 




/ \ -2e 5 p(>) 



(197) 



+j, (196) 

1 L 

2 ~ 2" ' 

"| + ^, (198) 
^7 

y-5L, (199) 
X C(2)-iL + V (201) 



-i 

e y ~> 2" ' 4' 
-l-^ + 2L-^, (202) 

^ + ^C(2)-7L + ^L 2 , (203) 

101 5C (2) + ^L-5L 2 , (204) 



6 " v ' 3 

3157 35 281 35 r2 , rtftr , 

^ + yC(2)--L + y L 2 , (205) 

-3C(2) ln(2) - lc(2) - §C(3) + L + ^C(2)L - ^L 2 

(206) 

5 + 12C(2) ln(2) + 11C(2) + yC(3) - 3L - 10C(2)L + ^L 2 

-jU 3 , (207) 

~ - 48C(2) ln(2) - Hc(2) - f C(3) + |^ + 36C(2)L 

-^L 2 + II 3 , (208) 

^ + 192C(2) ln(2) + Z^C(2) + 83C(3) - + 136C(2)L 

+ 113 l2 _ 43^ (2og) 
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a 2 E {0) 



156965 _ OA/rtNl , n . 12323^^ 
_-7B8C(2)h(2)-— «2). 

1387 _ 34 Q . 



629 >/oX 13319 



576 

+ 524 C (2)L-ifI^ + ^3. 



-C(3) + 



144 



(210) 




/'o 



-2e 



=-i j=i 



5 



where: 



(211) 



a 2 F M) 

a 2 F (_1) 
a r (2) 

a i* (3) 
a * (4) 

° F (5) 

a 2 F (0) 

a 2 F (0) 
a r (2) 

a 2 F (0) 
a _r (3) 

a 2 F (0) 



1 + -L, 

3 3 ' 
113 

--HL, 

809 170 

1 L, 

9 3 

14779 805 
~~36 3~ ' 

-I + I c(2) - l L ' 

» _ I C(2) _ l L _ 3 L ,, 
9 3 w 9 2 ' 

1123 



214 A/ . 574 r 149 r2 

-« 2) + - i -- i - 



?09397 + 11M 



144 



14345 
36 



1445 



£ 2 , 



(212) 
(213) 
(214) 
(215) 
(216) 
(217) 
(218) 
(219) 
(220) 
(221) 



where: 




-2e 



i=-2 j=l 



5 /nr(i) 

2/ J 



(222) 



a Ct (1) 
a ^ 

n 2 r (-2) 

a Cr (3) 

n 2 r (-2) 

a Lr (4) 



4 2 ' 

37 rr 

6 



(223) 
(224) 
(225) 
(226) 
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a Lr (5) - 


533 35 ^ 
24 2 ' 


(227) 


fl^G^ ^ = 


1 1 1 

_l_i C(2 )-iL--L 2 , 
4 sv ; 2 4 » 


(228) 


a 2 G<^ = 


3 + ic(2)-2L + iL 2 , 


(229) 


a%? = 


-f-fc(2) + 1 3L-^, 


(230) 


^ = 


517 „ . . . 19b r ^ r9 
— + 5C(2)-— L + 5L 2 , 


(231) 


a Ct (5) - 


19309 35 911 35 r9 
C(2)+ L L , 
48 2 s v ; 3 2 ' 


(232) 


a Ct (1) - 




(233) 


„2 r (0) _ 


-5 - 8C(2) - 7 -((3) -*L + U(2)L + \l 2 + ±L 3 , 


(234) 




^ 7 + 34<(2) + ™((3) + 39L 21((2)L ^L 2 L 3 , 


(235) 


„2 r (0) 


- 1 Jf- 4 >)-35C(3)-™ 1 i + 70 C (2)I + 3 «V 





n 2 r (0) 

a Lr (5) 



10 r3 



3397 3421 . . 245 >/ . 14389 r 

W + ^ C(2) + - C(3) + ^2- L 
6443 2 _ 35 s 

24 3 ' 



245C(2)L 



(236) 



(237) 



8 Expansion for Q 2 a 

We list, in this section, the expansion of the vertex diagrams around Q 2 = 0. 
Putting y = Q 2 /a, and keeping terms up to the order y 3 , we have: 



where: 




a 2 A { - 1} 
a 2 A { - l) 
a 2 A (_1) 



(h ■ k 2 ) 



l 

4' 



5 

72 



377 



21600 ' 



j = _i j =0 



/'0 



~ 2 £<,v, 

j"=o 



(238) 
(239) 



(240) 
(241) 
(242) 
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43 




1 



(282) 




240 

-1, 
2 



(283) 



(284) 



9 




211 



(285) 



3600 




1-|C(2), 



(286) 




19 




(287) 



216 




3613 




(288) 



108000 



9 Summary 

We have carried out a complete investigation of the scalar integrals associated to 
all the QED 2-loop vertex graphs, for on-shell electrons and arbitrary momentum 
transfer t = —Q 2 in the D-continuous regularization scheme. After identifying all 
the occurring Master Integrals (Mi's), we have written the linear, non- homogeneous 
differential equations in Q 2 satisfied by the Mi's, expanded them in e = (4 — D)/2 
and solved the equations by means of the method of the variation of the constants of 
Euler. The method requires the solution of the associated homogeneous equations. 
It turns out that all the associated homogenous equations are trivial, or became 
trivial after a suitable choice of the Mi's for the graph topologies involving more 
than a single MI; typically one had to solve a first order homogeneous differential 
equation with simple rational coefficients. 

The repeated integrations implied by Euler's method are immediately performed, 
in close analytic form, in terms of Harmonic Polylogarithms of increasing weight; 
the maximum weight occurring in the results presented in this paper was 4 (as in the 
case of the zeroth order term in e of the double cross topology). By further iterations 
of the approach, one could almost mechanically obtain any additional term in the e 
expansion of the Mi's. 

The explicit analytic evaluation of the QED vertex form factors in terms of the 
Mi's will be presented elsewhere. 
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A Propagators 

We list here the denominators of the integral expressions appeared in the paper. 



B 1-loop ingredients 

We recall in this appendix some useful results about 1-loop diagrams, fundamental 
ingredients for the 2-loop calculations, obtained with the method of differential 
equations. They appear in the 2-loop integrals which factorize in the products of 
two 1-loop integral; due to the possible presence of extra powers of 1/e in their 
coefficients, we give the results of the e expansion up to the second order in e. 

The case of the massive bubble is exaustively examined. The differential equation 
is presented and solved, as usual in the e — > expansion. 



V 2 

v 7 
v 8 
v 9 

V n 

£>13 
X>14 

T> 1B 

T>16 



k 2 

K 2 i 

(fa + fa) 2 , 

ipx-hf, 

(pa - h) 2 , 

[kl + a] , 

[k 2 2 + a] , 

[(ki + k 2 f + a] , 

[(pi ~ h) 2 + a] , 

[(P2 + h) 2 + a] , 

[(p 2 - hf + a] , 

[(Pi +P2- kif + a] , 

[(pi +P2 - + a] , 

[(Pi - h - k 2 ) 2 + a] , 

[(p 2 + h + k 2 f + a] , 

[(pi +P2 - fa - k 2 ) 2 + a] . 



(289) 
(290) 
(291) 
(292) 
(293) 
(294) 
(295) 
(296) 
(297) 
(298) 
(299) 
(300) 
(301) 
(302) 
(303) 
(304) 



B.0.1 Tadpole 
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where: 



A-i 
a 

Ao 
a 

Ax 
a 

M 

a 



i=-l 

-1, 
-1, 
-1, 
-1 . 



(305) 



(306) 
(307) 
(308) 
(309) 



B.0.2 Fully massive bubble 

The topology under consideration has one MI. We choose the scalar integral itself: 



F(e,a,Q 2 




{d D k} 



[fc2 + a ] [(Q - k) 2 + a] 



(310) 



The corresponding first-order linear differential equation is the following: 



dF(e,a,Q 5 



1 



2e] 



Q 2 (Q 2 + 4a)J 
(1 - e)\ 1 



2a 



F(e,a,Q 2 ) 
1 



[Q 2 (Q 2 + 4a) 



Tie, a] 



(311) 



where T(e, a) is the Tadpole. 

As in the cases previously discussed, we use our knowledge on the analytical 
behaviour of the solution in order to find the initial condition. In fact, Eq. ()311j) 
shows two possible singularities for the function F(e, a, Q 2 ), for Q 2 = and for Q 2 = 
—4a. Only the second, nevertheless, is indeed a singularity for F, corresponding to 
the physical threshold. Multiplying Eq. ()311|) for Q 2 and taking the limit Q 2 — > 0, 
we obtain: 

Fie, a, Q 2 = 0) = -±—^-T(e, a) . (312) 

a 

We look for a solution of Eq. (J311|) . with initial condition (J312|) . expanded in 
Laurent series around e = 0: 



F(e,a,Q 2 )= J] e'i^a, Q 2 ) + G(e 3 ) . 



(313) 



i=-l 



The associated homogeneous equation at e = is 



df(a,y) 
dy 



1 
2 



"1 1 

y (y + Aa) 



f(a,y) 



(314) 
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which has the following solution: 



f{e,a,y) = k x /l + ^-. 



(315) 



We can find the solution of the non- homogeneous equation, order by order in e, 
by means of the Euler's method of the variation of the constant k. We have: 



Fi{a,Q 2 



1 + 



4a 

~Q 2 



dy 



(y + 4a) 



1 + 



[A - A-i] 



1 n / \ 1/1 

_(y + 4a) 2a V?/ 



(316) 



where the coefficients Ai are those of Eqs. (|3U6H3UlJ|) . The determination of the 
constants of integration ki is made by imposing that the solution satisfies Eq. 
In terms of the variable x, defined in Eq. f!25|) . the solution reads: 




(4-D) 



{d D k} 



[k 2 + a] \{Q - k) 2 + a] 



V/W i= _ 1 



where: 

Bo 



Bo 



I- 

2-2 
4-4 



1 



1 



2 (1-x) 
1 1 



2 (1-x) 



#(0,x), 
C(2) 



o +H(0,x) + -H(p,0,x) 



-#(-1,0, x) 



4 + 4 



2 



1 + -#(0,x) -if(-l,x) 



(317) 

(318) 
(319) 



(320) 



C(3) 



H{0, x) + -H(0, 0, x) - H(-l, 0, x) + ^H{0, 0, 0, x) 



--i/(-l, 0, 0, x) - H(0, -1, 0, x) - (-1, -1, 0, x) 



(321) 



B.0.3 Bubble on the mass-shell 




/'o 



{4 - D) '{d^} 



fc 2 [(Q - k) 2 + a] 



Q 2 =-a 
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a(l - 2e) 



where: 



?=-i 

£-1 = 
i?o = 
^ = 

i?2 — 



1, 

2, 
4, 



(322) 



(323) 
(324) 
(325) 
(326) 



B.0.4 Scalar Vertex at 1 loop 




where: 
aF 



aFo 



1 



(27r)^- 2 ) A; 2 - A;) 2 + a] [(p 2 + k) 2 + a] 



(1 - 2 e ) 1 



e (Q 2 + 4a) 



(l-e) 1 



e (Q 2 + 4a) 



-e 2 



i=-l 



1 

2 



(1 + x) (1-x) 

1 



H(0,x), 



= -2 



(1 + x) (1-x). 
1 1 



{^-^ff(0,0,x) + ff(-l,0,x)|>, 



(1 + x) (1-x) 



C(3) C(2) 



-/f(0,x)-ff(-l,x) 



+\h(0, 0, 0, x) + ff (-1, -1, 0, x) - \h(-1, 0, 0, x) 



--ff(0,-l,0,x)|>, 



(1 + x) (1-x) 



9C 2 (2) 
20 



+ C(3)[ff(0,x)-2/f(-l,x)] 



+C(2) 



(327) 

(328) 
(329) 



(330) 



-if (0, 0, x) + 2/f (-1, -1, x) - H (-1, 0, x) - ff (0, -1, x) 
2 



— ff (0, 0, 0, 0, x) + 4/f (-1, -1,-1, 0, x) - 2fT(-l, -1, 0, 0, x) 
-2H(-1, 0, -1, 0, x) - 2H(0, -1, -1, 0, x) + 4/f (-1, 0, 0, 0, x) 
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+H(0, -1, 0, 0, x) + H (0, 0, -1, 0, x) 



(331) 



C Reducible 2-loop diagrams 

In this appendix we give the expressions of the reducible diagrams of Figs. (j3HHl)- 




2( 4 - D ) / SAD 



{d u h}{d D k 2 } 



1 



ViV 2 V 9 V 10 Vi 4 



(^E^+^e), 



where: 



(i) 



aE 



aE { }\ 



aE, 



(i) 



(1-x) (1 + x) 
1 1 



2 |_(1 -x) (1+x) 
+2#(-l,0,x)] 



H(0,x), 
[((2)-2H(0,x)-H(0,0,x) 

[2C(2) + 2C(3)-(4 + 3C(2))F(0,x) 



(332) 
(333) 

(334) 

(335) 



_(l-x) (1+x). 
-2C(2)if(-l, x) - 2H(0, 0, x) + 4ff(-l, 0, x) - 5F(0, 0, 0, x) 
+2H(0, -1, 0, x) + 2#(-l, 0, 0, x) - 4ff(-l, -1, 0, x)] . (336) 



where: 



(2) 




,, 2 ( 4 - D ) / SAD 



{d^fcl}{d D fc 2 } 



1 



V 4 V 6 V 7 V 14 V 15 



4) V + <^) 



(337) 
(338) 



(1-x) (1 + x) 



C 2 (2) 



+ 2C(3)/T(0,x)-C(2)fI(0,0,x) 



1 

~2 

-2C(2)i?(l, 0, x) - 4ff(0, 0, -1, 0, x) + 2#(0, 0, 1, 0, x) 

-4# (0, 1, 0, 0, x) - 2H(1, 0, 0, 0, x) 1 . (339) 




/'o 



2(4 " D) I {d D h}{d D k 2 } 



V 2 V 4 V 6 V 8 V U 



4) V + <^) 
Mo 



(340) 
(341) 
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where: 



aE, 



(3) 



1 



27C 2 (2) 



+ C(3)#(0,x) + 3C(2)#(0,0,x) 



(l-x) (1 + x) J { 10 
-6C(2)//(0, -1, x) + 2#(0, 0, 0, 0, x) + 2H(0, 1, 0, 0, x) 

-2/7(0, -1,0, 0,x) 




= ^ (4 - D) / {d D h}{d D k 2 } — — 1— - 
" j2 z l Ef + O (e) , 



/'o 



i=-2 



where: 



(l-x) (1 + x) 
1 1 



#(0,x), 



(342) 

(343) 
(344) 

(345) 



(4) 



1 



(1 +x) 

1 



(1-x) (1 + x) 

1- 1 



(1 + x) 



[C(2)-#(0,0,x) + 2#(-l,0,x)] (346) 
[6C(2) + 2//(0,0,x)]- 



(1-x) 



[C(3) - 2tf(0, x) - C(2)H(-l,x) - if (0, 0, 0, x) 



(1 + x). 

-2if(-l, -1, 0, x) + 0, 0, x) + if(0, -1, 0, x)] . (347) 



/i? 4 " D) / {</%}{<*%} 



V 2 V 6 V 8 V 12 V m 



V, "Effete), 

H'O, 



(348) 
(349) 



where: 



a^ 5) = 



(l-x) 



(l-x) J 



[3C(3) + 4#(-l, 0, 0, x) - 4#(0, -1, 0, x) 



+2#(0,l,0,x) -2#(l,0,0,x)]. 



(350) 



<0 



/ {d D h}{d D k 2 } 



V 1 V 5 V 7 V 8 V 10 



-2e 

a \ -,(6) 



(351) 
(352) 
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where: 



aE® = 6C(2)ln2-^C(3). 



(353) 



where: 



(7) 



aE 



aE, 



(7) 




2(4-D) / fjD 



{d u k 1 }{d D k 2 } 



V 1 V 7 V 9 V 10 V 13 



1 



(1-x) (l + x)_ 

1 1 



(1-x) (1 + x) 



+2H(-l,0,x)} + 
1 [ 1 



i=-2 

H(0,x), 

[((2)-2H(0,x)-H(0,0,x) 
1 



(1-x) 
1 



1 



H(0,0,x), 



(354) 
(355) 

(356) 



(1-x). 

.(1-x) (1 + x)J [2C(2) + 2C(3)-(4- C (2M(0,x) 

-2C(2)#(-1, x) + 2H(0, 0, x) - 4if (-1, 0, x) - H(0, 0, 0, x) 
-4if(-l, -1, 0, x) + 2i7(0, -1, 0, x) + 2if(-l, 0, 0, x) 

[C(2)iJ(0, x) -2H(0, 0, x) + 3#(0, 0, 0, x) 



(357) 



(1 — x) |_ (1 — x 
-4#(-l, 0, 0, x) - 2H(0, -1, 0, x)] . 



(358) 



where: 



v 12 v 16 



(8) 



E 



-2e 



j=-2 



2 
5 



(1-x) 



#(0,x), 



^ 8) = 



(1-x) 

'jihf H ^ x) 



[C(2)-5//(0,x) + 2//(-l,0,x)] 



(359) 
(360) 

(361) 
(362) 

(363) 
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where: 



(12) 



E 



(l-x) (1 + x) 

1 1 



2 [(l-x) (1 + x) 
+2H(-l,0,x)\, 



E, 



(12) 



H(0,x), 
[((2)-H(0,x)-H(0,0,x) 

(l-x) (1 + x)J [C(2) + 2 C (3)-(l-C(2))if(0 > x) 

-((2)H(-1, x) - H(0, 0, x) + 4ff(-l, 0, x) - #(0, 0, 0, x) 
-4if (-1, -1, 0, x) + 2#(0, -1, 0, x) + 2H(-1, 0, 0, x)] . 



where: 



2(4-D) IsjDuxsjD 



{d v k 1 }{d u k 2 } 



V 1 V 7 V W V 13 



e;(13) 
^-2 

^(13) 



(13) 



1, 

4- 
12 + 



1 - 



1 - 



(l-x) 
2 



tf(0,x), 



(l-x)J 
+2#(-l,0,x)]. 



[C(2)-4tf(0,x)-tf(0,0,z) 



where: 



-COr = 4 (4 - D) [ {d D k l} {d D k 2 } 



v 2 v A v & v 11 



/'o 



-2e 



£^+0(6), 



=-2 



1, 

4, 
12 
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where: 



where: 



»T~ D) f {d D md D h}^^- (392) 
J2 ^ (15) + O (e) , (393) 

i=-2 



a 
/'o 



R (15) 



Ei\ 5) 

a 
a 



-2e 



-1 , (394) 
-3 , (395) 
-7. (396) 



»T~ D) J{d D ki}{d D k 2 }^^- (397) 
£ ^if 6) +0(e), (398) 



-2e 

a 



/'o 



i=-2 



a 

£(16) 

a 



-1, (399) 
-3 , (400) 
-7. (401) 
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